Conformists and contrarians in a Kuramoto model with identical natural frequencies 
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We consider a variant of the Kuramoto model, in which all the oscillators are now assumed 
to have the same natural frequency, but some of them are negatively coupled to the mean field. 
These "contrarian" oscillators tend to align in antiphase with the mean field, whereas the positively 
coupled "conformist" oscillators favor an in-phase relationship. The interplay between these effects 
can lead to rich dynamics. In addition to a splitting of the population into two diametrically opposed 
factions, the system can also display traveling waves, complete incoherence, and a blurred version of 
the two-faction state. Exact solutions for these states and their bifurcations are obtained by means 
of the Watanabe-Strogatz transformation and the Ott-Antonsen ansatz. Curiously, this system 
of oscillators with identical frequencies turns out to exhibit more complicated dynamics than its 
counterpart with heterogeneous natural frequencies. 

PACS numbers: 05.45.Xt, 89.75.-k 



I. INTRODUCTION 

The Kuramoto model [T] of coupled oscillators has 
been used to shed light on many diverse systems of phys- 
ical interest, particularly those involving synchronization 
transitions. Examples include Josephson junction ar- 
rays [2], charge-density waves [3J, laser arrays @], col- 
lective atomic recoil lasers [5] , bubbly fluids [5] , neutrino 
flavor oscillations [7], electrochemical oscillators [3], and 
human crowd behavior [9]. 

Originally, however, the Kuramoto model had no 
known physical applications; these were only discovered 
years later. Kuramoto was led to his model solely by con- 
siderations of mathematical tractability. He was seeking 
an exactly solvable many-oscillator system displaying a 
phase transition to mutual synchronization, in hopes of 
illuminating this novel critical phenomenon seen earlier 
by Winfree in his simulations of biological rhythms |10j . 

In that same spirit, we have begun investigating a fam- 
ily of simple models that generalize the Kuramoto model 
in one key respect: they include both positive and neg- 
ative coupling in the same system. Positive coupling, 
analogous to a ferromagnetic interaction, tends to align 
the oscillators in phase. Negative coupling, analogous to 
an antiferromagnetic interaction, drives oscillators apart 
and favors a phase difference of ir. When both types of 
coupling are present, the system becomes frustrated. In 
this case very little is known about what sorts of dynam- 
ics and equilibrium states might follow. 

Even the mean-field version of such systems remains 
mysterious. Twenty years ago, in pioneering work, Daido 
found evidence that Kuramoto models with mixed pos- 
itive and negative coupling could undergo a glass tran- 
sition jllj . but the existence and properties of such an 
"oscillator glass" remain unclear [T^] . Other models with 
mixed attractive/repulsive interactions have since been 
explored by several authors, who were also motivated by 
analogies to spin glasses, as well as to neural networks 



with mixed excitatory and inhibitory connections [13] . 
In each instance it has been difficult to understand the 
behavior of these models because of their inherent nonlin- 
earity, quenched random interactions, and large numbers 
of degrees of freedom. 

Inspired by Kuramoto's success in explaining Win- 
free's synchronization transition through the use of a 
toy model, we wondered whether Daido's oscillator glass 
transition might be similarly rationalized by studying 
much simpler models with mixed coupling. In this paper 
we analyze the behavior of one such model and find, un- 
fortunately, that this particular simplification does not 
exhibit an oscillator glass. Nevertheless, this negative 
result still provides valuable information. It shows that 
certain types of frustration are insufficient to produce 
an oscillator glass, and thereby constrains the possible 
mechanisms at work. 

Furthermore, the model does display some interesting 
new dynamical phenomena, as we discuss below. And 
although we are unaware of any physical realization of 
the model studied here, we suspect that such realizations 
may exist, given the model's simplicity, and given the 
history of the Kuramoto model itself, whose physical rel- 
evance was established only after the model had been 
proposed on theoretical grounds. 

The governing equations for the model are 
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k=l 

( s) 

where <f>j is the phase of the jth oscillator in the s- 
subpopulation, lo is its natural frequency, and N is the 
total number of oscillators. The oscillators in subpopu- 
lation 1 are assumed to have positive coupling (K i > 0) 
to all the other oscillators in the system, whereas those 
in subpopulation 2 have negative coupling [K% < 0). 

Equation ([I]) differs from the classic Kuramoto model 
in that the distributed natural frequencies Uj have been 
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replaced by a uniform natural frequency w, and the sin- 
gle positive coupling constant K has been replaced by 
a two- valued coupling constant K s . In an earlier paper, 
we considered the case in which Uj was kept heteroge- 
neous [T3]. As we will see here however, the long-time 
dynamics is actually more complicated for the homoge- 
neous case. This finding is consistent with previous stud- 
ies of identical oscillators (see Refs. (TSIIH] for example). 

What is unusual about this model is that its pairwise 
interactions need not be symmetric. For example, oscil- 
lator k could be coupled positively to oscillator j while j 
is coupled negatively to k in return. This leads to a novel 
type of frustration. Although unfamiliar, it may be phys- 
ically realizable in certain kinds of series arrays of Joseph- 
son junctions [2] or in liquid crystal spatial light mod- 
ulators suitably coupled by global optoelectronic feed- 
back [T7]- Because of its asymmetry, this form of cou- 
pling is non-variational; no energy function exists, and 
the dynamics do not correspond to relaxation or gradi- 
ent descent down an energy landscape. While this might 
seem unnatural in some physical settings (e.g., magnetic 
spin systems), it is more plausible in certain social or po- 
litical contexts. In particular, if we set uj = (as can 
be done without loss of generality by going into a suit- 
able rotating frame, or equivalently, by replacing (f>j with 
(f>j +u)t), the model starts to resemble some of the existing 
models of social opinion formation |18j . 

To see the connection, imagine a spectrum of opin- 
ions or attitudes that can be laid out as points on a cir- 
cle, rather than as points on a line. For instance, Bin- 
more [IH] has argued that political attitudes are more 
properly represented this way than as the usual linear 
continuum from left wing to right wing. Now consider 
a population of indifferent individuals who have no pre- 
ferred phase along the circle, or, to continue the analogy, 
no inherent political preference. All that matters to them 
is what other people think. Such an individual updates 
his or her political "phase" continuously, based on where 
he or she stands in relation to the prevailing sentiment. 
Some individuals - the conformists - want to be in phase 
with conventional wisdom, whatever it happens to be, 
whereas contrarians reflexively oppose it. 

The question is: depending on the relative proportions 
of conformists and contrarians, and depending on how 
intensely they react to the prevailing opinion, what will 
this population do in the long run? Split into two camps? 
Fail to reach any consensus at all? Or cycle through 
all attitudes periodically? As we'll see, all of these are 
possible long-term outcomes, depending on the choice of 
model parameters. 

Let p denote the fraction of oscillators with positive 
coupling; thus the system consists of pN conformists and 
qN contrarians, where q = 1 — p (Fig. [IJ. When p = 0, 
all the oscillators repel one another, a case explored in 
Ref. [20] • 

In what follows, we will examine the dynamics of this 
system as its parameters are varied. We will continue to 
use the metaphorical language of conformists and con- 




FIG. 1: A schematic plot of a set of eight identical oscillators 
with positive and negative coupling. Open circles represent 
oscillators with negative coupling to all the others; filled cir- 
cles denote oscillators with positive coupling. When p = 0, all 
oscillators are contrarians. However, as p is increased, con- 
formists begin to appear and eventually replace all contrarians 
as p approaches 1. 



trarians, although perhaps we should stress that we do 
not intend the model to be taken literally as a descrip- 
tion of real social situations. It is a toy model. Like the 
Kuramoto model itself, it is being offered on theoretical 
grounds, without any particular physical realization in 
mind. The goal is to clarify the dynamical consequences 
of mixed coupling, by investigating a particularly simple 
and tractable special case. Our hope is that such an in- 
vestigation may bring us a step closer toward solving the 
puzzle of "oscillator glass" [TTl - fT3] . 



II. DIMENSIONAL REDUCTION 

The dynamical system given by Eq. ([I]) enjoys a highly 
non-generic structure. It has N — 6 constants of motion, 
for all N > 6. In geometrical terms, its phase space is fo- 
liated by an (N — 6)-parameter family of six-dimensional 
invariant manifolds. 

These results can be seen explicitly by using a theoreti- 
cal device discovered by Watanabe and Strogatz [21] , and 
recently generalized by Pikovsky and Rosenblum |22) . 
Rewrite Eq. ([I]) as 



f + K s g cos 



»<•> 



where / = 



Ke h sin 



(1/N) X^feLi cos 0fe- As discussed above, assume that 
lu = without loss of generality. Watanabe and Strogatz 
showed that all solutions (f>j (i) of Eq. (|2| can be expressed 
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tan 
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where the ipi*' in Eq. &3) are constant, and j s (t), Vf^i), 
and Q s (t) evolve according to the ordinary differential 
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equations 

is = 

* s = 

Q s = 



-(1 - j*)K s (gsme a - hcose s ), 
— — -K s g cos O s + h sin O s 



(4) 



/ cos 8 S + h sin Q s 



Here again the variables (p^it) denote the oscillator 
phases in the s-subpopulation, where s = 1,2. The 

(s) 

constants %fn represent a set of fixed phases on which 
the transformation operates. For example, if we set 
7s (0) = * s (0) = 9 S (0) = 0, then the ipj are just the 
initial phases 4>j(0). Since the variables 7 s (i)> ^ s (t), and 
0s (i) are the same for all j within each subpopulation, 
the flow governed by Eq. Q is effectively 6-dimensional, 
as claimed. It describes the dynamics restricted to the 
invariant manifold labelled by the choice of the constants 
ipj , of which N — 6 turn out to be independent [21] ■ 

A further reduction is possible in the continuum limit 
N — > oo, in the special case where the phases ipj are 
uniformly distributed around the circle. Then the trans- 
formation Eq. ([3]) maps old phases ipj to new phases (pj 
such that for each subpopulation, a uniform distribu- 
tion of the ipj maps onto a Poisson kernel distribution 
of (pj [15, 21-23 . This implies that the set of states in 
which each subpopulation is distributed like a Poisson 
kernel is dynamically invariant (see Ref. [IH1 EH |2"TH2"3"] 
for more about this). From here on, we will refer to this 
distinguished invariant manifold as the Poisson subman- 
ifold. 

Incidentally, these considerations underlie the (other- 
wise seemingly miraculous) ansatz discovered by Ott and 
Antonsen [IB]. They found that Poisson kernels are also 
dynamically invariant for the original Kuramoto model, 
where the oscillator frequencies arc non-identical. This 
beautiful invariance property has its origin in group the- 
ory [23j and has allowed many new insights to be gained 
into the dynamics of the Kuramoto model and its rela- 
tives [El [24H26] . 

On the Poisson submanifold, two of the equations in 
the system Q decouple from the other four. Thus, as we 
will see in detail below, the dynamics become effectively 
four-dimensional there. And because of an additional 
rotational symmetry (stemming from the fact that the 
right hand side of Eq. ([I]) involves only phase differences, 
not absolute phases), the flow can be further reduced to 
a three-dimensional system, which appears later in this 
paper as Eq. 



III. SIMULATION OF THE REDUCED SYSTEM 

We now numerically explore the dynamics of the six- 
dimensional system given by Eq. Q to get a sense of 
its equilibrium states. To do so, we recall that p is the 



fraction of the N oscillators that are conformists, and we 
define C as the relative intensity of the conformist cou- 
pling: C = K\j{K\ — K 2 ). Thus, values of C close to 
1 mean the conformists are much more intense in their 
desire to be like everyone else, as compared to the rel- 
atively mild obstinacy of the contrarians. On the other 
hand, when C is close to 0, the conformists are tepid 
while the contrarians are passionate. 

For our initial conditions, we choose each 7 S (0), ^(O) 
and O s (0) uniformly at random from [0, 1), [— tt,tt) and 
[— 7r, 7r), respectively. In addition, we set the N constants 
tpj such that each subpopulation is evenly spaced on the 
interval [—an, got) for a < 1. For example if the pN con- 
formists are indexed first and the qN contrarians after, 
we set ipj to be 



4>j = 



2a7r(j - pN/2) 
p~N 

2an(j -pN - qN/2) 
qN 



j = 1, 



,pN, 



(5) 



j = P N + 1,...,N. 



As we noted in the previous section, choosing a = 1 con- 
fines the trajectories to a distinguished submanifold of 
the phase space in which the new phases <pj for each sub- 
population are distributed like a Poisson kernel. Choos- 
ing a < 1 instead gives an initial condition off this special 
manifold and therefore allows the system to explore other 
parts of phase space [3TJ [57] . 

We begin by numerically integrating Eq. Q from ini- 
tial conditions on the Poisson submanifold. After tran- 
sients have decayed we compute the final phase density 
P{<p) and order parameter Z , defined by 



Re i8 = 



1 N 

3=1 



(6) 



We also compute the final order parameters of each sub- 
population: 



_L v 



(7) 



where J\ = {1,. . . ,pN} and J 2 = {pN+1, . . . , N}. Here 
r s represents the degree of synchronization of the subpop- 
ulation s, and 9 S denotes its average phase. Likewise, N s 
is the number of oscillators in this subpopulation. The 
integration itself is done using Heun's method with a time 
step of 0.01. 

Out of the whole 10 5 time steps, the first 7 x 10 4 time 
steps are discarded as transients, after which the quan- 
tities of interest were measured and averaged for the re- 
maining time steps. The system generally seems to end 
up in one of four states (Fig. [2]) : 

(a) The incoherent state, in which both the conformists 
and contrarians are uniformly distributed around 
the unit circle in the complex plane, yielding r\ = 
?'2 = 0. In terms of the political analogy discussed 
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(a) incoherent state 



(b) blurred n-state 



P(+) 



P(+) 
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(c) traveling wave state 










(d) Tr-state 







FIG. 2: (Color online) Four states commonly observed in the 
long-time behavior of numerical trials on the Poisson subman- 
ifold. The integration was performed for TV = 10 4 oscillators 
and the final states are presented as histograms with a bin 
size of 0.01. Conformist oscillators are shown in blue and 
contrarian oscillators in red. The four states shown are (a) 
the incoherent state at (p,C) = (1/20,2/3), (b) a blurred 
vr-state at (p, C) = (1/4,2/3), (c) a traveling wave state at 
(p,C) = (1/4,2/3) and (d) a Tr-state at (p, G) = (4/5,2/3). 
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FIG. 3: Trajectories of the order parameter Z(t) both (a) on 
and (b)-(d) off the Poisson submanifold. The specific param- 
eter values (p,C,a) for these states are (a) (1/2,2/3,1), (b) 
(2/5,2/3,1/2), (c) (1/2,2/3,1/2) and (d) (11/20,2/3,3/10). 



earlier, this means that no predominant attitude 
emerges in the population. All points on the polit- 
ical spectrum are equally represented. 

(b) A one-parameter family of blurred ir-states, corre- 
sponding to non-uniformly distributed populations 
of conformists and contrarians on the unit circle. 
The peaks of their phase distributions are blurred 
and separated from one another by an angle of n. 
Here the political interpretation is that two main 
factions have emerged, in diametrical opposition to 
one another. They could lie anywhere on the polit- 
ical spectrum, but once they emerge, the contrar- 
ians oppose the conformist view. And because of 
the blurred nature of both peaks, "fringe" views 
are also present on either side of the two main at- 
titudes. 

(c) A traveling wave state, in which the conformists 
and contrarians exhibit full and partial phase syn- 
chrony, respectively, with the peaks of their phase 
distributions offset by an angle less than ir. Here 
the conformists are unified in their views, yet that 
consensus view keeps changing, periodically cycling 
through all possible points on the political spec- 
trum. Meanwhile the contrarians oppose them, but 
not quite diametrically, and their opinions remain 
dispersed throughout. 

(d) The TT-state, in which the conformists and contrari- 
ans are completely synchronized into two antipo- 
dal delta functions (and thus ri = r-i = 1 and 
|#2 — 0i\ = tt). This simple state represents im- 
placable polarization between two unified and un- 
changing factions. 



The offset by an angle less than tt for the traveling wave 
state induces a nonzero wave speed. Hence Z traces out 
a circular orbit, as shown in Fig. [3^ a). Interestingly, the 
traveling wave state has been also found in analogous 
systems with heterogeneous natural frequencies [14]. 

The long-time dynamics of Z becomes substantially 
more complicated when we evenly space the constants 
tpj on the interval [—air, air) for a < 1, corresponding 
to initial conditions lying off the Poisson submanifold. 
Typical trajectories appear either quasiperiodic or pos- 
sibly chaotic in these cases, as shown in Figs. [3](b)-(d) . 
Similarly non-periodic behavior off the Poisson subman- 
ifold has been seen in other systems of oscillators with 
identical frequencies [T51 HH H2] . 

IV. ANALYSIS OF THE REDUCED SYSTEM 

According to Refs. [El [27], the Watanabe-Strogatz 
transformation is the real part of the Mobius transfor- 
mation. Additionally, we can convert 7 S in Eq. Q to r s 
via the relation j s = — 2r s /(l + r 2 s ) to obtain 

f s = 1 -^K s Re{Ze~ m °), 

* s = ^±K a Im{Ze- ie °), (8) 

Q s = l -^K s lm{Ze- i@ '), 
zr s 

where Re and Im denote the real and imaginary parts of 
their arguments. We here note that the relation between 
7 S and r s works only on the Poisson submanifold. It is 
not satisfied elsewhere (see details in Ref. [H]). Using 
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the fact that Z — pZ\ + qZ 2 and defining <5 
Eq. Q becomes 



r-2 



C(l 

-(1 

sin 8 



r\){pr\ + qr 2 cos 8), 

- rl){pri cos 8 + qr 2 ), 



(9) 



p(l - G) — + rir 2 
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By a fixed point analysis of Eq. ([£]), we can show that 
the four states found above by simulation are in fact the 
only generic equilibrium states of the reduced system re- 
stricted to the Poisson submanifold [57]. We summarize 
several interesting points of this analysis in the remain- 
der of this section and compute the order parameter R 
for the four different states. We do the latter by mak- 
ing repeated use of the relation Z = pZ\ + qZ 2 . De- 
composed more fully, this is Z = pr\e lBl + qr 2 e 102 or 
Z = (pn + qr 2 e lS )e tf>1 . 



A. Incoherent state 

We start with the easiest case: The incoherent state 
has T\ = r 2 = 0, so its order parameter R is zero. By 
a linear stability analysis [37J, we find that this state is 
stable when p < 1 — C . This gives us our first bifurcation 
value of p: p^ = 1 — C. 



B. Blurred 7r-states 

The one-parameter family of blurred 7r-states is given 
by the following fixed points of Eq. ^9|: 8 — n and 
pr\ = qr 2 (where r±, r 2 7^ 0) . By our above equations 
for Z , this implies R — 0. Linear stability analysis then 
shows that the blurred 7r-states farthest from the incoher- 
ent state begin to lose stability at p a = (1 — y/2C — l)/2, 
whereas loss of stability for the entire set of blurred tt- 
states occurs as p nears pt = 1 — C. Hence, there are 
stable blurred 7r-states on the same region that the inco- 
herent state is stable. 



C. Traveling wave state 

Next we turn to the traveling wave state. From 
Fig. [2^c), it is clear that the conformists are fully syn- 
chronized (ri = 1) for this state, so Eq. (J9j) reduces to 



r-2 =-(l-C)(l-r|)(pcos<y + gr 2 ), 



sin 8 



p(l - C) 



(H 



- 2qCr 2 



(10) 



where sin<5 7^ 0. If we then solve for the fixed points of 
this system, we obtain: 



'2 



p(l - C) 



2qC-p(l-C) 



8 : 



^{-q^/p). (11) 



This solution can only exist for < r\ < 1, which implies 
that the traveling wave state exists only for p greater than 
p a = (1 — V2C — l)/2 and less than p c — C. Within this 
region of existence, we can determine the order parameter 
R for the traveling wave state using the fact that Z = 
prie 101 + qr 2 e l ° 2 . Since by definition Z — Re 10 , we know 
that R 2 = ZZ, or 



R 2 = p 2 + 2pqr 2 cos 8 + q 



2„2 



Substituting Eq. (11) into Eq. (12) gives 



R=\P 2 - 



p(l- P ) 2 (l-C) 



2C-p(l + C) ' 



(12) 



(13) 



which goes to zero at p — p a , as the numerical data 
suggests. 



D. The vr-state 

Lastly, we consider the 7r-state. As Fig. [2jd) suggests, 
the subpopulations are both synchronized (rx = r 2 = 1) 
and antipodal to each other (8 — 71"). Substituting r\ = 
r 2 = 1 and 8 = n into Z = pZ\ + qZ 2 yields 



R = 2p—1 



(14) 



which can only be positive for p > 1/2. A more system- 
atic stability analysis of the 7r-state shows that it is stable 
for p > max{C, 1/2} [57]. For example, when the con- 
formist coupling is twice the magnitude of the contrarian 
coupling (e.g. K\ = 1 and K 2 = —1/2, and so C = 2/3), 
Eq. ( 14 ) implies that the 7r-state is stable for p greater 
than p c = 2/3. Interestingly, this agrees with the p c be- 
yond which the traveling wave state does not exist. Yet 
p c is not a bifurcation point since 8 is always it for the 
7r-state but does not approach 7r for the traveling wave 
state as p approaches p c . Instead, the r% = 1 nullcline 
corresponding to f 2 = and a parallel nullcline corre- 
sponding to 8 = approach each other as p approaches 
p c until they coincide and form a line of fixed points at 
P = PcW\- 

The above analysis indicates that the stable states 
of the reduced system reach their boundaries of sta- 
bility (and sometimes also existence) at three transi- 
tional points: p a = (1 — \J2C — l)/2, p& = 1 — C, and 
Pc = max{C, 1/2}. We can verify this by first comput- 
ing R for the numerically discovered fixed points found 
at various values of p, and then plotting the theoretical 
curves that we found above on top of these numerical 
data. The results in Fig. [4] illustrate the level of agree- 
ment between simulation and theory. 

We finish by summarizing which states lie in which 
regions delimited by the transitional points p a , pi, and 
p c - In general, there are four regions of behavior, and for 
C = 2/3, their boundaries arep a = (l-y/l/3)/2 « 0.21, 
Pt, = 1/3 and p c — 2/3. On these regions, we have found 
the following states: 
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FIG. 4: (Color online) Behavior of the order parameter R 
for the four different types of equilibrium states as a function 
of the conformist fraction p. Parameter values: C — 2/3, 
TV = 10 4 oscillators. The colored symbols denote R values 
of fixed points found in numerical simulations of Eq. ([I]) , for 
initial conditions on the Poisson submanifold. Blue asterisks 
denote both the incoherent state and blurred 7r-states; open 
magenta squares denote the traveling wave state; filled cyan 
squares denote the 7r-state. The solid red curve traces the 
theoretical value of R for the traveling wave state where it 
exists, and the blue line gives the theoretical R for the n- 
state on its domain of existence. 

I: On < p < p a , both the incoherent state and all 
of the blurred 7r-states are stable. 

II: On p a < p < p , the incoherent state, some of the 
blurred 7r-states, and the traveling wave state are 
stable. 

Ill: On p < p < p c , only the traveling wave state is 
stable. 

IV: On p c < p < 1, only the 7r-state is stable. 

To be more precise, the notion of stability being used 
above is that of stability within the Poisson submanifold, 



not stability within the full phase space. All states within 
the Poisson submanifold are neutrally stable to pertur- 
bations off the submanifold, because such perturbations 
carry the system onto another invariant manifold of the 
foliation discussed earlier. 



V. SUMMARY 

In this paper, we considered a system of identical os- 
cillators with positive and negative global coupling and 
investigated how the interplay between the positive and 
negative interaction affected the collective dynamics and 
equilibrium states of the system. We reduced the dynam- 
ics of our system from N dimensions to six by means of 
the Watanabe-Strogatz transformation, and found that 
in the infinite- N limit there are four types of equilibrium 
states of the system on a special submanifold of the phase 
space (the invariant manifold of phase distributions given 
by Poisson kernels). Using both numerical and analyti- 
cal techniques, we characterized each of these equilibrium 
states, paying particular attention to the illustrative case 
in which the conformists were coupled twice as strongly 
to the mean field as the contrarians were. Even for this 
slice of parameters, however, we found that a menagerie 
of complicated states exist throughout the phase space. 



Acknowledgments 

H.H acknowledges the hospitality of Cornell Univer- 
sity. This research was supported by the Chonbuk Na- 
tional University (H. H.) and NSF Grants CCF-0835706 
and CCF-0832782 (S. H. S.). We thank Seth Marvel for 
several useful discussions and the Korea Institute for Ad- 
vanced Study for providing computing resources for the 
project. 



[1] Y. Kuramoto, Chemical Oscillations, Waves and Turbu- 
lence (Springer, New York, 1984). 

[2] J. W. Swift, S. H. Strogatz, and K. Wiesenfeld, Physica 
D 55, 239 (1992); K. Wiesenfeld, Physica B 222, 315 

(1996) ; K. Wiesenfeld, P. Colet, and S. H. Strogatz, Phys. 
Rev. Lett. 76, 404 (1996); ibid., Phys. Rev. E 57, 1563 
(1998); K. Park and M. Y. Choi, Phys. Rev. B 56, 387 

(1997) ; G. Filatrella, N. F. Pedersen, and K. Wiesenfeld, 
Phys. Rev. E 61, 2513 (2000); B. C. Daniels, S. T. M. 
Dissanayake, and B. R. Trees, Phys. Rev. E 67, 026216 
(2003). 

[3] G. Griiner and A. Zettl, Phys. Rep. 119, 117 (1985); G. 
Griiner, Rev. Mod. Phys. 60, 1129 (1988); S. H. Strogatz 
and R. E. Mirollo, J. Phys. A 21, L699 (1988); ibid., 



Physica D 31, 143 (1989); C. M. Marcus, S. H. Strogatz, 
and R. M. Westervelt, Phys. Rev. B 40, 5588 (1989). 

[4] Y. Braiman, T. A. B. Kennedy, K. Wiesenfeld, and A. 
Khibnik, Phys. Rev. A 52, 1500 (1995); S. Yu. Kourtcha- 
tov, V. V. Likhanskii, A. P. Napartovich, F. T. Arecchi, 
and A. Lapucci, Phys. Rev. A 52, 4089 (1995); G. Kozyr- 
eff, A. G. Vladimirov, and P. Mandel, Phys. Rev. Lett. 
85, 3809 (2000); ibid., Phys. Rev. E 64, 016613 (2001); 
R. A. Oliva and S. H. Strogatz, Int. J. Bifurcation Chaos 
11, 2359 (2001); A. G. Vladimirov, G. Kozyreff, and P. 
Mandel, Europhys. Lett. 61, 613 (2003). 

[5] C. von Cube et al, Phys. Rev. Lett. 93, 083601 (2004); J. 
Javaloyes, M. Perrin, A. Politi, Phys. Rev. E 78, 011108 
(2008). 



7 



[6] G. Russo and P. Smereka, SIAM. J. Appl. Math. 56, 327 
(1996). 

[7] J. Pantaleone, Phys. Rev. D 58, 073002 (1998). 
[8] N. Mazouz, K. Krischer, G. Flatgen, and G. Ertl, J. Phys. 
Chem. B 101, 2403 (1997); I. Z. Kiss, W. Wang, and J. 
L. Hudson, Chaos 12, 252 (2002); M. Wickramasinghe 
and I. Z. Kiss, Phys. Rev. E 83, 016210 (2011). 
[9] Z. Neda, E. Ravasz, Y. Brechet, T. Vicsek, and A.-L. 
Barabasi, Nature 403, 849 (2000); Z. Neda, E. Ravasz, 
T. Vicsek, Y. Brechet, and A.-L. Barabasi, Phys. Rev. E 
61, 6987 (2000). 

[10] A. T. Winfree, J. Theor. Biol. 16, 15 (1967); The Geom- 
etry of Biological Time (Springer, New York, 1980). 

[11] H. Daido, Phys. Rev. Lett. 68, 1073 (1992). 

[12] H. Daido, Phys. Rev. E 61, 2145 (2000); J. C. Stiller and 
G. Radons, Phys. Rev. E 58, 1789 (1998). 

[13] H. Daido, Prog. Theor. Phys. 77, 622 (1987); L. L. 
Bonilla, C. J. P. Vicente, and J. M. Rubi, J. Stat. Phys. 
70, 921 (1993); D. H. Zanette, Europhys. Lett. 72, 190 
(2005). 

[14] H. Hong and S. H. Strogatz, Phys. Rev. Lett. 106, 054102 
(2011). 

[15] S. A. Marvel and S. H. Strogatz, Chaos 19, 013132 
(2009). 

[16] E. Ott and T. M. Antonsen, Chaos 18, 037113 (2008); 
Chaos 19, 023117 (2009). 



[17] E. A. Rogers et al, Phys. Rev. Lett. 93, 084101 (2004). 
[18] P. L. Krapivsky and S. Redner, Phys. Rev. Lett. 90, 

238701 (2003); S. Galam, Physica A 333, 453 (2004); M. 

S. Lama, J. M. Lopez, and H. S. Wio, Europhys. Lett. 

72, 851 (2005). 

[19] K. Binmore, The breakdown of social contracts, in Social 
Dynamics, edited by S. N. Durlauf and H. Peyton Young, 
MIT Press, 2004. 

[20] L. S. Tsimring, N. F. Rulkov, M. L. Larsen, and M. Gab- 
bay, Phys. Rev. Lett. 95, 014101 (2005). 

[21] S. Watanabe and S. H. Strogatz, Physica D 74, 197 
(1994). 

[22] A. Pikovsky and M. Rosenblum, Phys. Rev. Lett. 101, 
264103 (2008). 

[23] S. A. Marvel, R. E. Mirollo, and S. H. Strogatz, Chaos 

19, 043104 (2009). 
[24] E. A. Martens, E. Barreto, S. H. Strogatz, E. Ott, P. So, 

and T. M. Antonsen, Phys. Rev. E 79, 026204 (2009). 
[25] L. M. Childs and S. H. Strogatz, Chaos 18, 043128 

(2008). 

[26] C. R. Laing, Chaos 19, 013113 (2009); Physica D 238, 
1569 (2009); Phys. Rev. E 81, 066221 (2010). 

[27] S. A. Marvel, H. Hong, and S. H. Strogatz (to be pub- 
lished) . 



